An Update on the Life Analysis 
of Spur Gears 

John J. Coy*, Dennis P. Townsendt, and Erwin V. Zaretskyt 

Gears used in aircraft and other applications may fail from scoring, tooth fracture due to 
bending fatigue, or surface pitting fatigue. Figure 1 illustrates a few of the failure modes that have 
occurred during experimental gear testing at the Lewis Research Center. The scoring type of failure is 
usually lubrication related and can be corrected by proper lubricant selection and/or changes in gear 
operating conditions (ref. 1). Tooth breakage is caused by tooth loads that produce bending stresses 
above the endurance limit of the material (ref. 2). It is usually accepted that the endurance limit, if it 
does exist, can be predicted from available stress-life (S-N) curves for the material being used (ref. 3). 
Current methods (ref. 4) of predicting gear surface-pitting failures are similar to those used for 
predicting the bending fatigue limit. According to the method of reference 4, the maximum surface 
contact stress (Hertz stress) should be limited to a value less than the surface endurance limit of the 
gear material. It is commonly believed that the gears would then have an infinite surface-pitting life. 
But, based on gearing studies (refs. 5 and 6) and on rolling-element bearing-life studies (ref. 7), there 
is no real evidence to support the concept of a surface-fatigue limit under normal operating 
conditions of bearings and gears. Rather, it appears that all gears, even if designed properly to avoid 
failure by scoring and tooth-bending fatigue, will eventually succumb to surface pitting in much the 
same way as rolling -element bearings. 

In references 8 to 10 a method of predicting the surface-fatigue lives of low-contact -ratio spur 
gears was presented. The method was based on the commonly accepted Lundberg-Palmgren theory 
that has been used for many years to predict the lives of rolling-element bearings. Fatigue testing of 
spur gears of various materials, using different methods of manufacture and lubricants, has been in 
progress for the past 10 years at Lewis. 

The objective of the information presented in this paper is to summarize and update the spur- 
gear life analysis and to relate the analysis to some of the experimental results from NASA tests. 


Symbols 

b half width of Hertzian contact, m (in.) 

C gear center distance, m (in.); constant of proportionality, see eq. (14) 
c orthogonal shear stress exponent 

D a rolling-element diameter, m (in.) 

E Young’s modulus, Pa (psi) 

e Weibull exponent 

/ face width of tooth in contact, m (in.) 

Gjo 10 percent life of a gear, millions of revolutions 
h depth to critical stress exponent 

K proportionality constant 

L\q 10 percent life, millions of revolutions (Mr), or, equivalently, hr 
Lj gear rotations in millions 

/ involute profile arc length, m (in.) 

M contact ratio 

m g gear ratio, N 1 /N 2 

‘Propulsion Laboratory, U.S. Army Research and Technology Laboratories (AVRADCOM), NASA Lewis Research Center. 
tNASA Lewis Research Center. 
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NASA TEST GEAR 



SCORING 



FRACTURE FROM 
SURFACE PIT 



FRACTURE AT 
TOOTH ROOT 



SURFACE PITTING 


Figure 1. - Typical gear tooth failure modes. From testsrun on NASA Lewis 
spur gear tester. 


N number of teeth 

Pb base pitch, m/tooth (in. /tooth) 

Q load normal to involute profile, N (lb) 

Q c dynamic capacity for bearing, N (lb) 

Q m dynamic capacity of gear mesh, N (lb) 
q maximum Hertz stress, Pa (psi) 

r pitch circle radius, m (in.) 

r a addendum circle radius, m (in.) 

r b base circle radius, m (in.) 

S probability of survival 

Tio 10 percent life of a tooth, millions of cycles 

U number of contact cycles per bearing revolution 

u stress cycles per revolution 

V volume, m3 (in. 3) 

x,y Cartesian coordinates, m (in.) 

Z contact path length, m (in.) 

z depth of occurrence of critical shearing stress, m (in.) 

0 roll angle increment, rad 

5 precontact roll angle, rad 

t) millions of stress cycles 

6 base circle roll angle, rad 

/x Poisson’s ratio 

p curvature radius, m (in.) 

Ep curvature sum, m - 1 (in. - *) 

r maximum subsurface orthogonal reversing shear stress, Pa (psi) 

tp pressure angle, rad 

Subscripts: 

H high load 

/ index representing 1 or 2 
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L low load 

M mesh of pinion and gear 

1 reference to high speed member, pinion 

2 reference to low speed member, gear 


Theory 


Theory of Gear Tooth Life 

The fatigue-life model proposed in 1947 by Lundberg (ref. 1 1) is the commonly accepted theory 
to determine the fatigue life of rolling-element bearings. The probability of survival as a function of 
stress cycles is expressed as 


108 s“ 


T c rj e V 
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Hence, if the probability of survival is specified, the life rj for the required reliability can be 
considered a function of the stressed volume V, the maximum critical stress t and the depth to the 
critical shearing stress z. As a result, the proportionality can be written as 
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The above formula reflects the idea that greater stress shortens life. The depth below the surface z at 
which the critical stress occurs is also a factor. A microcrack beginning at a point below the surface 
requires time to work its way to the surface. Therefore, we expect that life varies by an inverse power 
of depth to the critically stressed zone. 

The stressed volume Vis also an important factor. Pitting initiation occurs near any small stress- 
raising imperfection in the material. The larger the stressed volume, the greater the likelihood of 
fatigue failure. By the very nature of the fatigue-failure phenomenon, it is the repetitive application 
of stress that causes cumulative damage to the material. The greater the number of stress cycles, the 
greater the probability of failure. By experience it has been found that the failure distribution follows 
the Weibull model. 

As an expression for the stressed volume let 


VocfzJ 
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where / is the involute length across the heavy load zone as defined in the appendix. The choice of this 
length is a simplification, in that a more complicated treatment would be to integrate across the entire 
involute length, which is composed of both heavy and light load zones. The simplification is 
justifiable, however, since operation under lighter loads greatly diminishes the probability of failure. 
Another reason for choosing the heavily loaded zone is that all observed failures on the NASA gear 
tests occurred in that zone, mostly where the Hertz load on the pinion is greatest (ref. 12). It is also 
assumed that the most severe Hertz stress generated in the contact of the teeth is acting over the entire 
load zone of single-tooth contact. The calculation of stresses is considered next. From Hertz theory 
(ref. 13) the maximum stress at the surface where line contact is assumed is given by 



(4) 


The semiwidth of the contact strip is 


b = 


/ 8Q(1 -g2) 
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where Up is the curvature sum, such that 


Lp = — + — 


Pi P2 


( 6 ) 
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The expression for radius of curvature, p, is given in the appendix. According to reference 11, the 
cause of fatigue flaking of bearings is the maximum reversing orthogonal shear stress that occurs at a 
depth z below the surface and has an amplitude that varies between ±t where 


z = 


(7) 


and 


T = 


( 8 ) 


Using the foregoing expressions for stress, stressed volume, and depth of stress allows equation (1) to 
be written as 


log - a Q( c -h + i)/2f-(c- 
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Assuming a probability of survival of 90 percent and making use of the foregoing expressions for z, 
t, and V, we may express life as 


rl ocQ-(c-h + \)/2ef(c-h-l)/2e'Ep-(c + h-l)/2ei-l/e 
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Evaluation of Exponents 

It has been shown in reference 14 that life is inversely proportional to the 4.3 power of load and 
that the dispersion in life (Weibull exponent), e, equals 2.5 for AISI 9310 steel gears. From reference 
1 1 the dependence of dynamic capacity on bearing size leads to the following proportionality: 

Q c oc£>£' + ' , “ 3)/(c_/,+ 1) (11) 

where Q c is the dynamic capacity and is defined as the load that gives a bearing life of 1 million stress 
cycles at a 90-percent reliability level and where D a represents ball size and as such is a measure of the 
size effect. In reference 11 the exponent of D a is 1.07. Using the foregoing values, c and h are 
calculated to be 23.2525 and 2.7525, respectively. The life equation then becomes 

r? = ATQ-4.3/3.9i: iO -5/-0.4 (12) 

where A" is a constant of proportionality and t) is tooth life in millions of stress cycles. Based on the 
data presented in reference 14, it has been calculated that A' = 3.72x 10 18 when the units used in the 
equation are pounds and inches. 

Gear Life 

According to the relation for survival probability of a single gear tooth under constant service 
conditions, 

log ^ ocj f (13) 

For a 90-percent survival rate S = 0.90 and for the corresponding life v~Tio, a constant of 
proportionality C is defined by 

1 Og 0 ^ 9 =Cr > 0 (14) 


The life T \ o is conceptually identical to the 5 10 life for rolling-element bearings. For any general 
survival rate the following equation holds: 
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From basic probability theory the probability of survival for N teeth is the product of survival 
probabilities for the individual teeth. Then, the probability of survival for the gear with TV, teeth 


Sj = S N i 


io 4, =n -(^) 108 os 


where L, is gear rotations and u is stress cycles per gear revolution. Normally u = 1 , but in the case of 
a bull gear used to collect power from a number of inputs, u would be set equal to the number of 
inputs, assuming that equal power was applied through each input. For a 90-percent survival for the 
gear S, = 0.9, L, = Gi 0 , and 


108 55 


Then the 10-percent life for the gear is related to the 10 percent life of a single tooth by the following 
equation: 




For the case of an idler gear, where power is transferred without loss from an input gear to an output 
gear through the idler gear, both sides of the tooth are loaded once per revolution. The stressed 
volume is doubled in this case, which is the same as doubling the number of teeth. 




For most circumstances the basic tooth life, as calculated by equation (12), can be applied to 
teeth on both member gears in a mesh. However, for cases where the tooth proportions of one gear 
member are different from those of the mating member, a separate calculation for each gear tooth 
must be made. 


Life and Dynamic Capacity of the Mesh 

From probability theory the probability of survival of the two gears in mesh is given by 
S = SyS 2 (21) 

and, therefore, 

(22) 

where the subscript 1 denotes the higher speed gear (pinion) and 2 denotes the other gear. This 
equation is valid for a 90-percent probability of survival for the mesh combination of gear 1 and gear 
2. In this equation it should be noticed that all lives L should be expressed in the same time base, for 
example, in hours or in rotations of gear 1 . 
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The dynamic capacity Q m for the mesh is defined as the normal tooth load that may be carried 
with a 90-percent probability of survival for one million revolutions of the highest speed shaft. Then 
for any other pitch line load, the corresponding mesh life is calculated from the following equation: 


^io- 



(23) 


For example, imagine that gear 1 is on the higher speed gear shaft. Also assume that the teeth on the 
gear members are of the same proportions and are running in a simple mesh. In this case, based on 
the foregoing analysis, the dynamic capacity of the gear mesh is given by 

f[2.Sf9.15 ”1 0.093 

(24) 

Zpi2.5[N x {\ + m\-5) A 



Results and Discussion 

Gear fatigue tests were conducted at the Lewis Research Center with gears made from AISI 9310 
steel (ref. 14). 

The test rig is shown in figure 2. The dimensions of the gears are given in table I. 

Three groups of vacuum arc remelted (VAR) AISI 9310 case carburized and hardened steel gears 
were fatigue tested under loads of 463 x 1 0 3 , 578 x 10 3 , and 694 X 10 3 N/m (2645, 3305, and 3966 
lb/in.), which produced maximum Hertz stresses of 1.531, 1.710, and 1.875 GPa (222 000, 248 000, 
and 272 000 psi). The lubricant was a superrefined naphthenic mineral oil with a 5 percent extreme- 
pressure additive package. Failure of the gears was due to surface fatigue pitting. Test results were 
statistically evaluated using the methods of reference 15. The results of these tests are plotted on 
Weibull coordinates in figure 3. Weibull coordinates are the log-log of the reciprocal of the 
probability of survival graduated as the statistical percent of specimens failed (ordinate) against the 
log of stress cycles to failure or system life (abscissa). In each test group there were 19 failures out of 
19 tests. All failures occurred in the zone of single tooth loading at or just below the pitch diameter. 

The 90-percent confidence interval limits were determined for each group of test data (fig. 3). 
The interpretation of these limits is that the true life determined from an extremely large sample of 
gears running at each stress condition will fall between these confidence limits 90 percent of the time. 
Were these confidence limits to overlap, the life differences as determined from the test would not be 
considered statistically significant. The consistent trend of decreasing life with increasing stress 
indicates good statistical significance in the data that resulted from these three test series. 


TABLE 1. - SPUR GEAR BATA 


[Gear tolerance per ASMA class IE.] 


Tooth width in contact, cm (in.) 

Tooth width, cm (in. ) 

Number of teeth 

Diametral pitch 

Circular pitch, cm (in.) 

Whole depth, cm (in. ) 

Addendum, cm (in. ) 

Chordal tooth thickness reference, cm (i.n. 

Pressure angle, deg 

Pitch diameter, cm (in.) 

Outside diameter, cm (in.) 

Root f i 1 let , cm ( in. ) 

Measurement over pins, cm (in.) 

Pin diameter, cm (in.) 

Backlash reference, cm (in.) 

Tip rel ief , cm ( in. ) 


0.25 (0.10) 

0.64 (0.25) 

28 

8 

0.9975 (0.3927) 

0.762 (0.300) 

0.318 (0.125) 

) 0.485 (0.191) 

; 20 

8.890 (3.500) 

9.525 (3.750) 

. . 0.102 to 0.152 (0.04 to 0.06) 

9.603 to 9.630 (3.7807 to 3.7915) 

0.549 (0.216) 

0.0254 (0.010) 

0.001 to 0.0015 (0.0004 to 0.0006) 
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(b) SCHEMATIC DIAGRAM. 


Figure 2. - NASA Lewis Research Center's gear fatigue test apparatus. 


The American Gear Manufacturers Association (AGMA) has published two standards for tooth 
surface fatigue, AGMA 210.02 and AGMA 41 1.02 (refs. 4 amd 16). AGMA 210.02 provides for an 
endurance limit for surface fatigue below which it is implied that no failure should occur. In practice, 
there is a finite surface fatigue life at all loads. AGMA 411.02 recognizes this finite-life condition. 
Therefore, it does not contain an endurance limit in the load-life curve but does show a continuous 
decrease in life with increasing load. Both AGMA standards are illustrated in figure 4. The AGMA 
load-life curves shown are for a 99-percent probability of survival or the L \ life. The experimental L\, 
L\q, and £50 lives from the data of figure 3 are plotted for comparison. 










(c) TANGENTIAL LOAD W,/f 694X10 3 N/M (3966 Ib/in. ). 


Figure 3. - Comparison of life-prediction theory with experimental results for 
VAR AISI 9310 steel spur gears. Speed, 10 000 rpm. 


It is evident that the load-life relation used by AGMA is different from the experimental results 
reported herein. The difference between the AGMA life prediction and the experimental lives could 
be the result of differences in stressed volume. The AGMA standard does not consider the effects of 
stressed volume, which may be considerably different from that of the test gears used herein. The 
larger the volume of material stressed, the greater the probability of failure or the lower the life of a 
particular gearset. Therefore, changing the size or contact radius of a gearset, even though the same 
contact stress is maintained, would have an effect on gear life. 
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Figure 4. - Comparison of experimental life for VAR AISI 9310 spur gears with 
AGMA life prediction. Speed, 10 000 rpm; lubricant superrefined naphthenic 
mineral oil with additive package. 

It should be mentioned here that the Weibull slope e was assumed to be independent of the stress 
level in the original work by Lundburg and Palmgren. There is some evidence to suggest that the 
Weibull slope is a function of the applied stress. This was noticed by Schilke (ref. 6); and the tests 
conducted at Lewis (ref. 14) at various stress levels gave Weibull slopes from 1.9 to 2.9. The Weibull 
slope e = 2.5 that is suggested for use in this report is an average value that will serve in most 
applications where the stress is not unusually high or low. When life calculations are made using the 
method outlined in this paper, the predicted life can be considered a reasonably good engineering 
approximation to what may be expected in a practical gear application. However, the theoretical 
prediction does not consider material and processing factors such as material type, melting practice, 
or heat treatment, nor does it consider environmental factors such as lubrication and temperature. 
All these factors are known to be extremely important in their effect on rolling-element bearing life 
(ref. 17). There is no reason why the effects used to determine bearing life should be significantly 
different from those used to determine gear life. Figure 5 is a collection of life data based on gear 
tests conducted at the Lewis Research Center. The lives are compared with a base-line life of unity for 
AISI 9310 gears. The mode of failure was pitting. More test data obtained with gear specimens under 
various test conditions and different materials and lubricants continue to be required to establish 
and/or affirm the material and processing factors and the exponents c, h, and e for gears, as they are 
used in various applications. However, the general methods presented herein explain and support the 
use of the statistical methods to predict spur-gear fatigue life. Table II gives a sample calculation for 
a simple mesh of two gears. 


MATERIAL 



Figure 5. - Summary of gear fatigue lives based on NASA spur gear tests with 
various materials. 
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TABLE II. - SAMPLE CALCULATION 


Symbol 

Definition 

Equation 

Result 

Q 

Normal load 

Gi 

ven 



1610 N (363 

lb) 

<P 

Pressure angle 






20° 

Nl 

Number of teeth on gear 1 






28 

«2 

Number of teeth on gear 2 






28 

f 

Face width 





2.79 mm (0.11 

in . ) 

r l 

Pitch radius 





4.45 cm (1.75 

in. ) 

r 2 

Pitch radius 





4.45 cm (1.75 

in. ) 

r al 

Addendum radius 





4.76 cm (1.88 

in. ) 

r a2 

Addendum radius 





4.76 cm (1.88 

in. ) 

r t: 

Base radius 





4.17 cm (1.64 

in. ) 

r b2 

Base radius 

1 

r 



4.17 cm (1.64 

in. ) 

Z 

Length of contact path 

V r al ' r bl + V r a2 “ 

r b2 ' (r l 

+ r^, )s in <? 

1.63 cm (0.641 

in. ) 

P'- 

Base pitch 

2nr b i/Ni 




9.35 mm (0.368 

in. ) 

3LI 

Roll angle increment 

(Z - Pb ) / r bl 




0.166 

rad 

BhI 

Roll angle increment 

(2p D - Z)/r bl 




0.058 

rad 

6 1 

Pre contact roll angle 

(r l + r 2 )sin * ' ( r a2 - 

V' 2 

b2) 


0.169 

rad 



r bl 






*1 

Length of involute du'-ing 
single tooth contact 

r bl S Hl( S l + 6 L1 + 7 

\ 

B Hl y 



0.879 mm (0.035 

in. ) 

•n 

Radius of curvature 

r bl ( S 1 + Bli) 




1.39 cm (0.549 

in. ) 

P 2 

Radius of curvature 

(r i + rjjsin <p - pj 




1.65 cm (0.648 

in.) 

l 0 

Curvature sum 

( 1 / o 1 ) * (l/o 2 ) 




8.55 cm -3 (3.36 in 

-* 1 ) 

T 10 

Tooth life 

3.72xl0 18 Q" 4 - 3 f 3 - 9 z o 

T 5 ,- 

-0.4 


59.6 

Mr* 

G I0 

Gear life 

T 2 5T 1/2 - 5 




15.7 

Mr* 

L 10 

Mesh life 

T 2 r"|- 1 /2.5 

L 2(1,G 10> J 




11.9 

Mr* 


*Mr = million revolutions. 


Summary of Results 

An analytical method for predicting surface fatigue life of gears was presented. General 
statistical methods were outlined, showing the application of the general methods to a simple gear 
mesh. Experimentally determined values for constants in the life equation were given. Comparison of 
the life theory with NASA test results and AGMA standards was made. Gear geometry pertinent to 
life calculations was reviewed. The following results were obtained: 
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1. The life analysis resulted in the following fundamental equation for gear tooth life: 


ij = ^Q-4.3/3.9E p -5/-0.4 


where K= 3.72 x 108 for calculations in pounds and inches. 

2. For a simple mesh of one gear with another the dynamic capacity is 


Qm = 


[ 


K2. 5/9.75 


Epl2.5/7Vj(l+ W l-5) 


] 0.093 


(12) 


(24) 


Appendix — Gear Geometry 

In figure 6 an involute tooth is shown. The involute curve may be thought of as the path traced 
by point A on the tangent line as it is rolled upon the base circle with radius r b . The radius of 
curvature of the involute at point A is given by p. 


p = r b d 

The x and y coordinates of point A are given by the following equations: 
x = /^(sin 6 - 6 cos 6) 
y = r b (cos 6 + 6 sin 6- 1) 

The differential element of involute arc length dl is 
dl = (dx 2 + dy 2) 1 72 = r b 6 d6 


(Al) 


(A2) 


(A3) 


The length of the path of contact Z is the distance between where the two addendum circles cross 
the line of action (as shown in fig. 7): 


A =\l r2 al- r ll + yJ r2 a2- r h- Cs[n * 


(A4) 



Figure 6. - Involute profile geometry. 
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where C is the distance between gear centers. The base pitch p b is defined as the distance from a point 
on one tooth to the corresponding point on the next tooth measured along the base circle (see fig. 7). 
It may also be thought of as the distance from one tooth to the next measured along the line of 
action. 


_ _ 2 * r bl 

Pb ~~N^~ 


2 *r b 2 
A 2 


(A5) 


The contact ratio M is defined as the average number of gear teeth in contact at one time. It is 
calculated by 

Z 

M- — (A6) 

Pb 


When two pairs of teeth are in contact, the load is assumed to be equally shared. For low- 
contact-ratio gears (1 <M<2 ) there are three load zones. The load zones for gear 1 are shown in 
figure 8. Two pairs of teeth are in contact during roll-angle increments & LU and only one pair during 
roll-angle increment j3#i. For gear 1 the roll-angle increment for which two tooth pairs are in contact 
is calculated by 

0Ll = (A7) 

r b\ 

And the roll angle movement for the heavily loaded zone is given by 

&H\ = (A8) 

r b\ 

Finally b\, the precontact roll angle is defined as the roll angle from the base of the involute to the 
start of the zone of action as shown in figures 8 and 9. The angle 6) is given by the following 
equation: 


6 ,= 


C si n <P~(rl 2 -rl 2 ) l/2 
?b\ 



CS-81-2493 


(A9) 


Figure 7. - Spur gears in mesh. 
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After integrating equation (A3) between the limits that bracket the heavy load zone, the result is the 
length of involute in the heavy load zone. 


h-fb\ 


+ Pli + 


) 


(A10) 


NORMAL 

LOAD 




Figure 3. - Load sharing diagram. Load on tooth 
for low-contact-ratio gear depends on roll 
angle. 


Figure 9. - Geometry for calculation of pre- 
contact roll angle, 6^. 
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